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Abstract 

Generation of high-energy proton bunch from interaction of an intense short circularly polar- 
ized(CP) laser pulse with a gas-hlled cone target(GCT) is investigated using two-dimensional 
particle-in-cell simulation. The GCT target consists of a hollow cone filled with near-critical gas- 
plasma and a thin foil attached to the tip of the cone. It is observed that as the laser pulse 
propagates in the gas-plasma, the nonlinear focusing will result in an enhancement of the laser 
pulse intensity. It is shown that a large number of energetic electrons are generated from the 
gas-plasma and accelerated by the self-focused laser pulse. The energetic electrons then transports 
through the foil, forming a backside sheath field which is stronger than that produced by a simple 
planar target. A quasi-monoenergetic proton beam with maximum energy of 181 MeV is produced 
from this GCT target irradiated by a CP laser pulse at an intensity of 2.6 x 10 20 W/cm 2 , which is 
nearly three times higher compared to simple planar target (67MeV). 
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I. INTRODUCTION 



With the rapid development of the chirped pulse amplification technique, generation of 
energetic ion beam by interactions of an ultra intense laser pulse with a solid target has 
become realizable. Such energetic ions can be promising for many scientific or societal 
applications, such as proton radiographyjlj], fast ignition for inertial confined fusion{2-4], or 
hadron-therapy|5]. For most of these applications, ion beams with high energy, low energy 
spread and high collimation are required. 



Depending on the target paraments and laser intensity, ions can be accelerated by several 
different mechanisms, such as shock acceleration j^, Q], light-pressure acceleration jg, 10-12], 



Coulomb explosion 



131 ]. target-normal sheath acceleration (TNSA) 



14-161, etc., as well as 



their combinations. In TNSA, the energetic electrons produced at the front of a target by 
the laser ponderomotive force propagate through the target into the backside vacuum can 
generate a sheath electrostatic field. The sheath field, of order 10 12 V/m, can accelerate the 
ions on the target back surface to high energies. However, the proton beams obtained in this 
way are typically characterized by low particle density, large divergence, and almost 100% 
energy spread. An improved TNSA scheme, using a microstructured double-layer (DL) 
target, can decrease the energy spread. The possibility to generate 1.3 MeV proton beams 
with energy dispersion~ 25% and 3 MeV carbon beams with energy dispersion ~ 17% using 
a microstructured DL target has already been demonstrated experimentally by Schwoerer 
et al.[l3] and Hegelich et al.jla]. respectively. 

A tiny hollow metal cone was first introduced in fast ignition experiments to shield the 
igniting laser pulse from the underdense region of the precompressed fuel plasma Ji| , and 
a remarkable increase in the thermal fusion-neutron yield was observed. Since then the 
cone target was intensively examined both in experiments and simulations [20IJ32 ] . PIC 
simulations showed that a cone target could nonlinearly guide and focus a laser beam, and 
improve the efficiency of the coupling and transport of the energy into dense plasma 33|, [34]. 
Accelerating proton beams using a cone target with open tip was also studied by Cao et 
al. (35!, and energetic ion bunches of high density were observed. In this paper, we report that 
quasi-monoenergetic proton beam with peak energy of 130MeV and maximum energy of 181 
MeV can be generated from a gas-filled cone target irradiated by a CP Gaussian laser pulse 
at an intensity of 2.6 x 10 20 W/cm 2 . The gas-filled cone target, as shown in Figfj], consists 
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FIG. 1: (color online). Schematic view of the interaction of a laser pulse with a GCT target. The 
left trapezoid represents the laser pulse (red), the trapezoid irradiated by the laser pulse represents 
the gas-plasma (blue), the two parallelograms besides the gas-plasma represent the cone (orange), 
the rectangle on the tip of the cone represents the foil (dark). 

of a hollow cone filled with near-critical heavy-ion gas-plasma and a thin foil attached to 
the tip of the cone. Our results indicate that, comparing with that from a simple proton 
target, energetic protons with smaller energy spread and higher energy can be obtained. 
This result can be attributed to the much higher electron density and temperature behind 
the foil and the small transverse size of the foil. The energetic electrons are generated from 
the gas-plasma and accelerated by the enhanced laser pulse, which undergoes self-focusing 
in gas-plasma and is even focused by the tip of the cone. These energetic electrons can 
easily propagate through the thin foil to form a stronger sheath field behind the foil (than 
that behind a planar target). Since the foil has a small transverse size, the protons in 
the foil are accelerated in the homogenous sheath field, so that the protons are accelerated 
longitudinally forward with smaller energy spread. In contrary to TNSA acceleration, here 
energetic electrons originate mainly from the gas rather than the solid target. 
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II. SIMULATION PARAMETERS 



We carried out simulation using a fully relativistic particle- in-cell code (KLAP2D) [8|,|9j.In 
simulations, the simulation box is 80 A x 20 A, where A = lfim is the laser wavelength, and 
contains 3200 x 800 cells. A CP laser pulse with a peak laser intensity of II = 2.6 x 
l0 20 W/cm 2 is normally incident from the left side, The pulse has a Gaussian radial profile 
with 2a = 10A full width at half maximum and a trapezoidal shape longitudinally with 
40T flat top and IT ramps on both sides, where T is the laser period. The corresponding 
peak dimensionless laser amplitude a = eE / (m e uc) is 9.8, where E, u, c,m e , and e are 
the laser electric field, frequency, speed of light in the vacuum, electron mass, and charge, 
respectively. The GCT target, as shown in FigJTJ consists of electrons, protons, and heavy 
carbon ions. The initial temperature of electrons, protons, and carbon ions is 10 eV. The 
cone has a width of lfim, and is located in 10 < z[fim] < 45 with the diameters of the left 
and right cone openings of 16/im and 2/im, respectively. For simplicity, the cone consists of 
carbon plasma with an electron density n e = 10n c , where n c = irm e c 2 / (e\) 2 is the critical 
density. The carbon gas-plasma is full in the cone with density n e = 0.8n c . The foil with 
2/im wide and 0.35/im thick is placed at z = 45/xm. It consists of a proton-carbon mixed 
plasma with an electron density n e = 40n c , and the ratio of C:H=1:1. 



III. SIMULATION RESULTS 



i a frequency u p smaller than the 



40] as soon as its total power P 



A laser beam propagating in underdense plasma wit 
laser frequency u undergoes relativistic self-focusing 
exceeds the critical value 

P cr « 17(lu/u p ) 2 GW; (1) 

The self-focusing is due to the relativistic mass increase of plasma electrons and the 
ponderomotive expulsion of electrons from the pulse region. Both effects lead to a local 
decrease of plasma frequency and an increase in refractive index. The strong non-linear 
self-focusing of the laser pulse propagating in the near-critical gas-plasma at t — 56T is 
shown in Figj2]^a). For clarity, only a part of the simulation box is shown. The spot size of 
the laser is focused to be smallest at z — 35/xm, and the transverse electric field is enhanced 
to 17 there, which is 1.7 times higher than the initial laser electric field. The pulse retains 
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FIG. 2: (color online). Left: Transverse electric field (in units of m e coj/e) at (a)i = 56T, (c) 
t = 64T. Right: (b) electron density (in units of critical density n c ), and (d) the longitudinal 



electron current J ze (in units of en c c) at t = 44T. 



its Gaussian radial profile, however, its spot size varies with the distance of propagation in 
a periodic manner. The smallest spot size at z — 35fim is about 3/im, while it varies to 
about Afim at z — 40fim. This result is due to dynamic balance between diffraction and 



non-linear self- 
approximation 



ocusing, which is also in good agreement with the analysis of the paraxial ray 



41] . When the laser propagates to the tip of the cone, it is even focused or 



squeezed by the tip of the cone for the small radius there, as is shown in FigfS^b). The spot 
size is focused to about 1/im at the tip of the cone (z = 45/im), with the transverse electric 
field as high as 20 there. The electrons that are initially at the front of the pulse are more 
efficiently accelerated as the pulse undergoes intensity enhancement due to self-focusing. 
Strong flows of relativistic electrons, axially comoving with the laser pulse, are observed 
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FIG. 3: (color online). (a) the proton energy spectrum of the protons behind the targets for different 
cases at t=120T.and (b) Evolution of the maximum proton energy for different cases. 

in the simulation, as shown in Figf2]^b)and Figfj^d). The maximum electron density near 
the axis is as high as 17n c , and the longitudinal electron current J Ze = —en e v Ze is about 
—I4en c c (negative J ze due to negative electron charge). These energetic electrons then 
transport through the thin foil and form a strong backside sheath field there. 

Figj3](a) shows the energy spectrum of the proton bunches behind the targets in the two 
cases at t — 120T. For the case shown, the maximum energy for the GCT target is about 
181 MeV, which is nearly three times higher than that of the planar target (65MeV) under 
the same conditions. The energy conversion efficiencies from laser to protons are 2.5% and 
0.7% for GCT and planar targets at t = 120T, respectively. For the GCT target, due 
to the small transverse size of the foil where the sheath field is homogenous, the energy 
spectrum of the proton bunch has a quasi-monoenergetic peak with energy dispersion of 
about 36%. In contrast, the energy spectrum from the planar target is much broadened due 
to multidimensional effects such as hole boring and other instabilities. The evolution of the 
maximum proton energy is shown in Figj3]^b). For the planar target, as the laser impinges 
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FIG. 4: (color online). (a)Evolution of the electrostatic fields at the place of the proton layer for the 
two cases of planar and GCT target and (b)the electron energy spectrum of the protons behind 
the targets at t=50T for planar target and t=80T for GCT target. 

on the target at t — 10T(the planar target is initially located at z — 10/xm), the maximum 
proton energy increases earlier than the GCT target (for which the laser impinges on the foil 
at t — 60T). However, as the electrostatic field at the place of the proton layer is much 
weaker for the planar target (see in FigJU^a)), the increase of proton energy is much slower 
than the GCT target. For the GCT target, the maximum proton energy increases rapidly 
from 3.3MeV to 120MeV in only 20T(from t = 60T to t = 80T), which is attributed to the 
strong electrostatic field during that time. At later time, the maximum proton energies in 
both cases remain almost constant. 

Evolutions of the electrostatic fields at the place of the proton layer for the planar and 
GCT targets are shown in FigJU^a), which explains the energy enhancement of GCT target 
in Figj3]^b). The electrostatic fields straight up quickly at t = 60T for the GCT target when 
the laser impinges on the target, because the energetic electrons generated from the gas- 
plasma reach the back side of the foil at t — 60T and establish a strong sheath field there(see 
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in FigfS^a)). The maximum electrostatic field is about 20 for GCT target at t = 60T, which 
is about 3.3 times higher than the planar target(6 at t — 44T). Since the energetic electrons 
expand away quickly, the electric fields decrease quickly after reaching the maximum for 
both cases. The electron energy spectrums behind the targets at t=50T for planar target 
and t=80T for GCT target are shown in FigJUb), at both times when the maximum energy 
of the electrons behind the targets is highest . We can see that the electron temperature 
and density are higher for the GCT target , which will result in a higher longitudinal field 
and eventually higher proton energy, as shown in FigJU^a) and Figj3]^a). 




t(T) t(T) 



FIG. 5: (color online). (a) Time evolution of the number of electrons behind the planar and GCT 
targets. (b)Time evolution of the number of electrons behind the planar originated from different 
places. 

Since the TNSA mechanism depends strongly on the charge seperation field established 
by the energetic electrons, it is of interest to investigate the electron number behind the 
target. FigJS^a) shows the time evolution of electron number behind the planar and GCT 
targets. For both cases, the electron numbers initially increase after the laser irradiates on 
the targets, and then flatten out. For GCT target, the electron number is nearly 5 times 
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higher than the planar target at t = 120T. From Fig|5^b) we can see the energetic electrons 
are almost generated from the gas-plasma (97% of the total number), while only a small 
number of the electrons are from the cone (3% of the total number). This result indicates 
that with the GCT target the efficiency of proton acceleration is determined by the electrons 
generated from the gas, which is quite different from the planar target (the electrons are from 
the target itself). 




FIG. 6: (color online). (a) Maximum proton energy for different gas-plasma density, (b) Maximum 
proton energy for different foil thickness. 

The effects of the gas-plasma density and the foil thickness of the GCT target are shown in 
FigJHJ It is found that the maximum proton energy remains almost the same(near 180MeV) 
while the gas-plasma density is between 0.4n c and 0.8n c , and the foil thickness is between 
0.3/im and 0.5/xm. These simulation results demonstrate that our acceleration scheme is 
robust. On the other hand, for the gas-plasma density, over-high gas-plasma density will 
result in much depletion of laser pulse, while over-low gas-plasma density will leads fewer 
enegertic electrons behind the foil, both will result in a decrease of the maximum proton 
energy. For the foil thickness, thin foils proved to be more efficient for ion acceleration in 
TNSA by hot-Electron recirculation j^], but over-thin foil will result in a quick expanding 
of the electrons, which will also result in a lower proton energy. 
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We have also stimulated the interaction of a Linear polarized(LP) laser pulse with a 
GCT target at the same intensity, while the other parameters are the same as in FigJTJ Our 
simulations results verify that similar phenomenon can also be observed with the LP laser 
pulse. A quasi-monoenergetic proton beam with peak energy of 139 MeV and maximum 
energy of 185 MeV can be generated. This indicates that our acceleration scheme can be 
also efficient for the LP laser pulse. 



IV. CONCLUSION 



In conclusion, proton acceleration from a GCT target is proposed to enhance the ion 
energy. A quasi-monoenergetic proton bunch with peak energy of 130MeV and maximum 
energy of 181MeV is achieved by using the GCT target at laser intensity of 2.6 x 10 20 W/cm 2 . 
It is nearly three times higher than that from the planar target. This result is attributed 
to a stronger electrostatic field behind the foil, which is formed by the energetic electrons 
generated and accelerated by the enhanced laser pulse in the gas-plasma. The effects of the 
gas-plasma density and the foil thickness have been investigated. The results demonstrate 
that our acceleration scheme is robust. Such GCT target may be difficult to make at present, 
however, with the rapid advance in nanofabrication technology such a small conical channel 



filled with gas-plasma should be realizable 



43| . Accordingly, the GCT target can remarkably 



reduce the cost of a laser driven ion accelerator in the applications such as cancer therapy. 
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